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Abstract 



> : Physic, aspects of a JINR project to reach .he planrred 5A GeV energy for the Au arrd U 

' beams and to increase the bombarding energy up to lOA GeV are discussed. The project aims 

■ to search for a possible formation of a strongly interacting mixed quark-hadron phase. The 
, relevant problems are exemplified. A need for scanning heavy-ion interactions in bombarding 

■ energy, collision centrality and isospin asymmetry is emphasized 
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^ I. INTRODUCTION 

o ' 

^ ' Over the last 25 years a lot of efforts have been made to search for new states of 

C ■ 

strongly interacting matter under extreme conditions of high temperature and/or baryon 
^ ! density, as predicted by Quantum Chromodynamics (QCD). These states are relevant 
' to understanding the evolution of the early Universe after Big Bang, the formation of 
neutron stars, and the physics of heavy-ion collisions. The latter is of great importance 
since it opens a way to reproduce these extreme conditions in the Earth laboratory. This 
explains a permanent trend of leading world research centers to construct new heavy ion 
accelerators for even higher colliding energy. 

Looking at the list of heavy-ion accelerators one can see that after the pioneering ex- 
periments at the Dubna Synchrophasotron, heavy-ion physics developed successfully at 
Bevalac (Berkley) with the bombarding energy to Eiab ~ 2A GeV, AGS (Brookhaven) 
Eiab ~ llA GeV, and SPS (CERN) Eiab ~ 160A GeV. The first two machines are 
closed now. The nuclear physics programs at SPS as well as at SIS (GSI, Darmstadt, 
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Eiab ~ lA GeV) are practically completed. The new relativistic heavy-ion collider (RHIC, 
Brookhaven) is intensively working in the ultrarelativistic energy range ^/snn ~ 200 GeV 
to search for signals of the quark-gluon plasma formation. In this respect, many hopes 
are related to the Large Hadron Collider (LHC, CERN) which will start to operate in the 
TeV region in two-three years. The low-energy scanning program at SPS (NA49 Collab- 
oration) revealed an interesting structure in the energy dependence of some observables 
at Eiab ~ 20 — 30A GeV which can be associated with the exit of an excited system 
into a deconfinement state. This fact essentially stimulates a new large project FAIR 
GSI (Darmstadt) for studying compressed baryonic matter in a large energy range of 
Eiab = 10 — 30A GeV which should come into operation after 2015 year [1']. These prob- 
lems are so attractive that the RHIC scientific society discusses a possibility to decrease 
the collider energy till the FAIR range [2]. 

On the other hand, in JINR there is a modern superconducting accelerator, Nuclotron, 
which has not realized its planned parameters yet. The Veksler and Baldin Laboratory of 
High Energy has certain experimental facilities and large experience in working with heavy 
ions. This study may actively be supported by theoretical investigations of the Bogoliubov 
Laboratory of Theoretical Physics. In ^ a program was proposed for investigating the 
dense strongly interacting QCD matter, formed in relativistic heavy ion collisions, based 
on acceleration of heavy ions like Au at the Nuclotron up to the maximal planned energy 
Eiab = 5A GeV. In view of new opened opportunities of the Nuclotron update to increase 
the bombarding energy up to lOA GeV and to get both Au and U ions with relativistic 
energies, the relevant physics problems are discussed in this paper. 



II. PHASE DIAGRAMS 



A convenient way to present a variety of possible states of strongly interacting matter is 
a phase diagram in terms of temperature T and baryon chemical potential fiB (or baryon 
density pb), as presented in Fig HI This picture shows in which region of the diagram the 
given phase is realized and which colliding energies are needed to populate this region. 

As is seen, a system, formed in a high energy collision, is fast heated and compressed 
and then starts to expand slowly reaching the freeze-out point which defines observable 
hadron quantities. At the Nuclotron energy Eiab = 5A GeV the system "looks" into 
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FIG. 1: Dynamical trajectories for central (6 = 2 fm) Au+Au collisions in T — (left panel) 
and T — fiB (right panel) plane for various bombarding energies calculated within the relativistic 
3-fluid hydrodynamics calculated with hadronic EoS 4]. Numbers near the trajectories are the 



evolution time moment. Phase boundaries are estimated in a two-phase bag model [&\. In the 
right panel the critical end-point calculated in the lattice QCD 6| is marked by the star and 
the shaded region corresponds to uncertainties of the bag model. 

the mixed phase for a short time (the left part of FigHj), however, uncertainties of these 
calculations are still large. To get agreement with the lattice data for finite T and /i^, 
masses of m, d quarks should be rather heavy, and the phase boundary is shifted towards 
higher hb 5| (the right part in Fig{T]). On the other hand, one should exercise caution 
because these dynamical trajectories were calculated for a pure hadronic gas equation of 
state, and the presence of a phase transition may noticeably change them. In addition, 
near the phase transition the strongly interacting QCD system behaves like a liquid rather 
than a gas, as was clarified recently at small fiB from both quark j3] and hadronic 8] 
side. As to high fiB values, it is a completely open question. One should also note 
that as follows from lattice QCD calculations for ^ the deconfinement temperature 
practically coincides with the transition temperature for chiral symmetry restoration while 
for baryon-rich matter it is still an open question. 

In FigJJl dynamical trajectories in the e — pb plane for the top project Nuclotron 
energy are given for two equations of state, without and with the first-order phase transi- 
tion 0,^]. The shaded band corresponds to the mixed phase of more restrictive EoS with 
heavy masses of u, d quarks, as shown in the right panel of Fig{Tl Nevertheless, both the 
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FIG. 2: Dynamical trajectory in the e — /fs^-plane for central Au+Au collisions calculated with 
two equations of state: pure hadronic (solid line) and with first order phase transition (dashed). 
Spatial averaging is done over the cube with the 4 fm sides and Lorentz contracted in the 
longitudinal direction ^. The shaded regions correspond to the mixed phase (upper one) and 
the non-reachable domain with the boundary condition T = 0, respectively. 

trajectories, being close to each other, spend some time in the mixed phase. The main 
difference between the results presented in FigJT] and Figl2] comes from the different aver- 
aging procedures used: in the first case, thermodynamic quantities were averaged over the 
whole volume of the interacting system, while in the second case, it was carried out only 
over a cube of 4 fm sides placed at the origin being Lorentz contracted along the colliding 
axis. Thus, for central collisions at the 5A GeV Nuclotron energy even if an average state 
of the whole strongly interacting system does not approach the mixed phase, an essential 
part of the system volume will spend a certain time in this mixed phase. An experimental 
consequence is that an expected observable signal of reaching the mixed phase should be 
rather weak. Note that for Eiab =10A GeV these conditions for a phase transition are 
fulfilled appreciably better. 

One should stress that the presented above dynamical trajectories and boundaries for 
the first-order phase transition were estimated for a system conserving a single charge, 
namely the baryonic charge. However, the behavior of this system near the phase transi- 
tion and particularly within the mixed phase will be qualitatively different if conservation 
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FIG. 3: Projections of the boundary surface on the {pB,T) plane at different x. The isospin ratios 
are x = 0, —0.1, —0.2, —0.3, —0.4, —0.5, starting from the right. The boundaries of the mixed 
phase with the quark and hadron phase are plotted as solid and dashed lines, respectively [l^ . 



of more than one charge is taken into account [10|, lUl]. So we turn to consideration of 
effects of additional conservation of the electric charge, or isospin of the system. 

We shall characterize the charge asymmetry by the electric-to-baryonic charge ratio 
Z/A = Pq/pb or by the isospin ratio x = {pq—pb)/2pb which are related as a; = Z/A—0.5. 
For Au+Au collisions we have Z/A = 0.4 or x = —0.1. 

An essential difference in the first-order phase transition for systems with one and two 
conserved charges stems from the fact that the phase boundary, which is a line for a single 
charge conservation, is getting a two-dimensional surface. In the latter case the conserved 
charges can now be shared by two phases in equilibrium in different concentration in each 
phase but consistent with the global charge conservation. If this is energetically favored 
by the internal forces and Fermi energies, then these degrees of freedom will be exploited 
by the system and will influence thermodynamic quantities. 

The boundary surface, so called binodal, may be parameterized in a different way, say 
as {T,pb,Pq}, {T,Pb,x}, {T,pb,x} or jp^pR xj. In Figll, some {T,pb,x} projections 



of a hadron-quark phase transition are shown 12|. 

It is seen that, for example, in symmetric matter at T = the baryon density ranges 
along 3.5po between the onset and completion of the transition. For x = the hadron 
boundary is close to that in the left panel of Fig{Tl This mixed phase domain becomes 
even larger in density for iso-asymmetric systems. Another important observation is that 
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the density at the onset (i.e. hadronic side of the phase boundary at the transition density 
Pc) decreases with increasing isospin asymmetry. If one compares points for x = and 
X = —0.2 at T = 0, a decrease is 6pB ~ 0.5po- But this effect is practically absent for 
T > 120 MeV. 

A more realistic description of the hadronic phase, which takes into account the den- 



sity dependence of hadron masses and coupling constants [IJ], was used for the results 
presented in FigHl The general trend of curves is quite similar to that in Figj3] but now 
the values of transition densities are hi ghe r though the same bag constant B^^^ =187 MeV 



was used in both the calculations 
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14l |. The main difference comes from the very low 



mass used for in-medium nucleon at the normal nuclear density and T = 0, = O.6M7V. 
In this case the coefficient C4 ahead the highest sigma-meson interaction term C^a"^ turns 
out to be negative and, therefore, such a system should be unstable iJ]. If the bag con- 
stant increases, the transition boundary moves to higher values. Note that in both the 
cases the temperature T{fiB = 0) is lower than the value 170-180 MeV expected from the 
lattice calculations. However, the lattice QCD data show that we deal with the first-order 
phase transition only for the high baryon chemical potential, above that for the critical 
end-point (see also the right panel in Fig{T]). Thus, it is not completely clear whether 
two-phase calculations should describe the lattice value of T(/iB = 0). Appropriate lattice 
QCD data are available only for small fiB and there are no data taking into consideration 
the conservation of both baryon and electric charges. 

Even noticeably stronger reduction of the transition density with increasing neutron 
fraction is predicted by the Catania group [l^ for T = 0. It was demonstrated that pc 
depends appreciably on poorly known properties of EoS at high baryon densities. This 
statement is illustrated in FiglSl Note quite low values of the transition densities which 
originate from small values of the bag constant used, i?^/^ = 140 — 170 MeV. This choice 
is argued by the following Witten hypothesis a state made of an approximately equal 
number of u, d, s quarks can have the energy per baryon number smaller than that for Fe; 
therefore, the quark matter is absolutely stable and purely quark stars may exist. This 
hypothesis put strong constraints on quark model parameters: The bag model parameter 
B should be small [18|. Today the existence of pure quark stars is not excluded but is 
rather considered as an exotic case. In this respect, the (non)observation of a strong 
isospin dependence of the transition density may be treated as a test of the Witten 
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FIG. 4: Hadron boundary of the mixed phase at different x [l^]- The isospin ratios are x = 
0, —0.1, —0.2, starting from the right. Hadronic phase is described within the relativistic mean- 
field approach with density-dependent hadron masses and couphng constants [ij]. The results 



are given for two values of the bag constant: B^^^ =187 (dashed lines) and 225 MeV (solid 
lines). 

hypothesis. 

The most striking feature of all results is a sharp decrease of the transition density 
Pc which takes place in the range Z/A ~ 0.3H-0.45, though the size of the reduction 
effect and its position on the Z/A axis, as is seen in FigJSl are strongly model dependent. 
Nevertheless, application of neutron-rich heavy ions seems to be very perspective to study 
the mixed QCD matter at energies even lower than the project Nuclotron energy 5A GeV. 
Dynamical trajectories presented in Figl5]show that the phase boundary may be reached 
at energies as low as lA GeV. 

As follows from FigslSHSl the most favorable temperatures for this reduction effect are 
in the range T < 80 MeV. It is well known that in central collisions of relativistic heavy 
ions the growth of the density ps > 3po with the energy increase is accompanied by the 
appropriate rise in temperature T > 100 MeV. Available stable nuclei cover only a very 
narrow region in isospin asymmetry Z/A ^ 0.39^0.40 exhausted by the ^^2^ and ^jlAu 
isotopes, respectively, what embarrasses checking the boundary reduction effect. The use 
of the long-lived ^^gAu isotope (r; ~150 days) allows one to move towards neutron-poor 
side of the boundary till Z/A =0.41. In addition, the construction of accelerators of intense 
beams of neutron-rich heavy ions with Z/A beyond this narrow region is an extremely 
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FIG. 5: Variation of the transition density with the proton fraction at T = for various EoS 
parameterizations of the relativistic mean-field theory: Dotted, dashed and solid lines correspond 
to the GM3 version to an additional inclusion of non-linearit y in p interaction [IT] and to 
an extra term describing the interaction of the isovector 5-meson 



15|, respectively. The points 

represent the dynamical trajectory in interaction zone during semi-central ^^^Sn + ^^^^n 



collisions at lA GeV (circles) and at 300A MeV (crosses) 
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complex technological problem. [15|. Some possibility to get lower temperatures and to 
extend the reached isospin asymmetry region is opened by study of semi-central rather 
than central collisions. This possibility is illustrated in Figj6] 

It is noteworthy that after about 10 fm/c the quadruple momentum is almost vanishing 
and a nice local equilibration is achieved. At this beam energy the maximum density 
coincides with reaching the thermalization. Then the system is quickly cooling while 
expanding. In Figj6]one can see that rather exotic nuclear matter is formed in a transient 
time of the order of 10 fm/c, having the baryon density around 3po; the temperature 
50 — 60 MeV , the energy density 500 MeV im~^ and the proton fraction between 0.35 
and 0.40. So the local neutron excess well inside the estimated mixed phase region may 
be even higher than that in colliding nuclei (note that for ^38^ ^-j^g isospin symmetry ratio 
is Z/A =0.387). 

The {p, pb, x} projections of EoS corresponding to the Maxwell construction are given 
in FigJT] for two isotherms. In accordance with the result familiar from the behavior of 
the systems with one conserved charge, for symmetric matter {x = 0) the pressure stays 
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FIG. 6: Time evolution of thermodynamic quantities inside a cubic cell of 2.5 fm wide, located in 
the center of mass of the system, is shown for semi-central 238^_238^ collisions at Eiab =1A GeV 
with b =7 fm. Baryon density, temperature, energy density and proton fraction are presented. 
Different curves in the upper-left panel are: black dots ~ the baryon density in po units: grey 
dots - the quadruple momentum in momentum space; squares - the resonance density jl5l ]. 



constant within the mixed phase. In contrast, at a; < the pressure changes during the 
transition increasing with the baryon density. This change of the pressure throughout the 
phase separation in asymmetric systems is an indication of a smoother transition than 
in symmetric systems as was first noted in the context of neutron star calculations 
Certainly, such behavior of pressure may crucially affect the evolution of two colliding 
nuclei. 

Along with pressure, the temperature as well as the baryonic and electric chemical 
potentials do not remain constant within the mixed phase. This behavior is responsible 
for disappearance of the entropy discontinuity at the give n p and x as pictured in the 
right panel of FigJHl This fact gave grounds for author 12] to claim that in the iso- 
asymmetric matter the first-order phase transition is smoothed and becomes the second- 
order phase transition. One should note that in spite of such behavior of the entropy, 
the first derivative of the thermodynamic potential with respect to temperature at the 



10 



400 ■ 



300 



200 



100 




300 



J— 
I 



100 



T 



/ 



0.5 1.0 



FIG. 7: Isoterms for different values of x at T = (left panel) and T = 50 
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121. For the = the latent heat is shown. 



constant chemical potentials (i.e. locally in the [li space) suffers a jump (see Fig|9l), what 
is evidence of the first-order phase transition in iso-asymmetric matter. At the same time 
if the {p, representation of entropy is used, the result of [I4] is well reproduced 

(see the left panel of Fig 18]). This mirrors a general property inherent to systems with the 
first-order phase transition: In the mixed phase the first derivatives of the thermodynamic 
potential have discontinuity in the {T, /ii, ...,/i„} representation but they are continuous 
in the {T, pi , p„} one. That is the direct consequence of the Gibbs phase-equilibrium 
conditions [10] from which follows that in the first representation the every point of 
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FIG. 9: Entropy per baryon as a function of temperature at the constant pressure for various 
isospin ratios (left panel) and entropy as a function of temperature at the constant chemical 
potentials fiB =1300, =300 and /ig =-100 MeV (right panel). The relativistic mean-field 
model with the density-dependent interaction is used for the hadronic phase and the bag model 
describes the quark phase [3]. 

the ra-dimensional mixed phase surface maps unambiguously on the {T, /xi, ...,/i„} space. 
Evidently, it is not the case in the {T, pi, p„} representation [19] as illustrated in 
FiglTOl Due to that the volume fraction of the second phase A = ^ jumps when the 
system enters in the mixed phase, if the first representation is used, while it does not 
in the second representation. Any extensive thermodynamic quantity A will follow this 
behavior of A since A = XAu + (1 — X)Aj. 

As was noted above, the presence of two conserving charges in the {T, I^q] repre- 
sentation changes dimensionality of the two-phase coexistence surface for the first-order 
phase transition. Namely, the one- dimensional line for the case of a single conserving 
charge transforms into a two-dimensional (or n-dimensional) (hyper) surface if two (or n) 
charges are conserved 1^, see FigJTOl Therefore, the manifold of critical points, that by 
definition is the boundary of the mixed phase where the system suffers the second order 
phase transition, also changes sufficiently: From two isolated points defining the limited 
line of the mixed phase it transforms into a one-dimensional curve {{n — l)-dimensional 
(hyper) surface] being a boundary of two-dimensional (ra-dimensional) coexistence (hy- 
per) surface |l9|. The different topologies of the mixed phase states may result in different 
important consequences. In particular, in the discussed now new projects of FAIR |l| 
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FIG. 10: Schematics view on the mixed phase boundaries in the {T, fi} (on left) and {T, p} 
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charges (middle) as well as the volume fraction (bottom) are presented 



and RHIC [2] aimed to search for manifestation of the critical end-point, being the point 
where the first-order phase transition of the QCD matter ends, it may turn out that it is 
not a point but hypersurface those dimensionality is defined by a number and properties 
of conserving charges |l9|. It is noteworthy that the location of the critical end-point 
in the T — fis plane was estimated in both various model and lattice QCD calculations, 
however, without taking into account isospin asymmetry degree of freedom. But even in 
this case the predictions vary wildly as demonstrated in FigJTT] taken from 20 1. 

The outstanding problems involved in understanding the hadron-quark matter phase 
transition, aside from the description of these phases themselves, are related to the geo- 
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FIG. 11: Theoretical (model and lattice) predictions for the location of the critical end-point. 
Two lines are obtained by lattice Taylor expansion where the lower curve corresponds to smaller 
quark masses (compare with Fig(T]). Points are calculated in different models specified in 
Errors/uncertainities are not shown. 

metric structure of the mixed phase and its evolution with varying the relative fraction A 
of phases. If one of the conserved charges is the electric charge, a geometric structure in 
the mixed phase is expected. The nucleation mechanism for cluster formation (say, quark 
drops in the hadron environment) is dominant in the metastable region of the first-order 
phase transition and formation of the mixed phase influenced by finite size effects due to 
nonvanishing surface tension at the interface between the hadronic and quark matter and 
the Coulomb energy of the formed drops. The geometric structure of the mixed phase was 



studied in some detail for neutron star matter (see the review-article [2l|]) but for nuclear 
matter uncertainties are rather large. In FigJT2l an estimate of the finite-size effects is 
presented for the nuclear case jisl]. 

As is clearly seen, the value of the density at which the mixed phase can be reached 
becomes larger, if the finite size effects are taken into consideration. This increase due to 
finite-size effects is larger for smaller values of the bag constant B (compare with results 
in FigHj). Nevertheless, a strong dependence of the transition energy on Z/A survives 
even in this case Il5l. 
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FIG. 12: Transition density as a function of ZjA at T = 50 MeV [15[1. The solid lines are 
obtained neglecting finite-size effects. The results plotted as long and short-dashed lines take 
into account these effects and correspond to B^l'^ = 160 MeV and B^l'^ = 170 MeV, respectively. 



Besides the isospin there is yet another nuclear parameter which may influence evo 



lution of a colliding system. It is the nuclear shape. It was noted earlier [22|, |23| that 
deformation and orientation affect compression, elliptic flow and particle production for 
collisions of Uranium nuclei. As seen from FigJT3| the compression in the tip-tip U+U 
collisions is about 30% higher and the region with p > 5po lasts approximately 40% longer 
than in the body-body collisions or spherical U+U collisions. Moreover, the nucleon el- 
liptic flow has some unique features which in principle may allow one to disentan gle t hese 
|. Such situation is valid for the energy range 1-20A GeV |22(] but 



two orientations 



uncertainties are still large |23| 



III. ON SIGNALS AND PRECURSORS 



Similarly to the general situation with searching for a quark-gluon plasma in ultrarel- 
ativistic collisions, there is no single crucial experiment which unambiguously solves the 
problem. It is quite evident that direct information on the mixed phase may be obtained 
only by means of weakly interacting photon and lepton probes. Sensitivity of global 
hadron observables to possible phase transitions is expected to be weak, but it might not 
be the case for more delicate characteristics. In any case, due to the proximity of the 
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FIG. 13: Evolution of central baryon density for Au+Au and U+U collisions at 20A GeV (left 
panel) and impact parameter dependence of nucleon elliptic flow at lOA GeV (rig ht panel) for 



different orientations of deformed Uranium nuclei as well as for spherical ones 
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phase diagram region under discussion to the confinement transition and chiral symmetry 
restoration, some precursory phenomena cannot be excluded at the bombarding energy 
below lOA GeV, which opens a new perspective for physical investigations at the Dubna 
Nuclotron. 

Properties of hadrons are expected to change in hot and/or dense baryon matter 
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25| . This change concerns hadronic masses and widths, first of all those for the a-meson 
as the chiral partner of pions which characterizes a degree of chiral symmetry violation 
and can serve as a "signal" of its restoration as well as the mixed phase formation. Rare 
decays in matter of vector mesons (particularly p and uj) are also very attractive. 

Nevertheless, there are theoretical proposals to probe chiral symmetry restoration in 
the vicinity of the phase transition boundary. In particular, it was shown H] that a 
two-photon decay of the cx-meson formed as an intermediate state in vrvr scattering may be 
a very attractive signal. As depicted in FigHU at temperature in the vicinity of the phase 
transition, when mo- ~ 2m,r, there is an anomalous peak in invariant mass spectra of 7 
pairs which may serve as a signal of the phase transition and formation of a mixed phase. 
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FIG. 14: Invariant mass spectra of 27 at = 80 MeV and different temperatures 
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FIG. 15: Invariant mass distribution of pairs of 7 quanta with the energy > 100 MeV after 
subtraction of the combinatorial background in the d+C 7+7+X reaction at the bombarding 
energy 2A GeV l28|. 

Certainly, there is a huge combinatorial background due to 7r° — > 77 decays, but the 
Nuclotron energy is expected to have some advantage against higher energy accelerators 
because the contribution of gamma's from deconfined quarks-gluons will be negligible. 

To measure this signal, a photon spectrometer is needed. The installation of this kind, 
two- arm 7-spectrometer PHOTON-2, is available at JINR. The ability of this installation 
in discriminating two correlated photons is demonstrated in FigJTSl 

The presented invariant mass spectrum of 77 correlations is based on 1.5 ■ 10^ selected 
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dC-interactions. A peak corresponding to about 5000 r^-meson decays is clearly seen at 



28( 1 ■ The measured rj mass 



Mr, = 540.5 ± 2.1 MeV with the resolution a = 33.4 ± 6.0 
is consistent with the table value. A huge peak from the vr^ 77 decay is strongly 
suppressed here and is not visible in the figure due to the selection criteria of events to 
increase the signal/background ratio. It is of interest that some additional structure with 
a maximum in the range of 340-460 MeV is seen which was not observed in analogous 

n 

experiments with lower statistics [28]. The nature of this resonance structure is under 
discussion now. 

Fluctuations in particle number are inherent in the first-order phase transition and for- 
mation of a metastable mixed phase. However, to a certain extent, peculiarities of a phase 
transition will be washed out by transition dynamics and subsequent hadron interactions 
which drive the system closer to equilibrium, which means the loss of information on the 
system prehistory. Nevertheless, it was understood that the study of fiuctuations in rela- 
tivistic strongly interacting matter may help with solving the problems mentioned above. 
Experimental data on event-by-event fluctuations (e.g., fluctuations in particle multiplic- 
ity, electric, baryon and strangeness charges) in nuclear collisions give a unique possib ility 



to test recent ideas on the origin of fluctuations in relativistic interacting systems |29l. |30|. 
Among them the suppression of event-by-event fluctuations of electric charge was pre- 
dicted [29!] as a consequence of deconfinement. Theoretical estimates of the magnitude 
of the charge fiuctuations indicate that they are much smaller in a quark-gluon plasma 
than in a hadron gas. Thus, naively, a decrease in the fiuctuations is expected when 
the collision energy crosses the threshold for the deconfinement phase transition. How- 
ever, this prediction is derived under assumptions that initial fiuctuations survive through 
hadronization and that their relaxation times in hadronic matter are significantly longer 



than in the hadronic stage of the collision 



29 



ail. 



More realistic estimates have been done recently for the hadron resonance-gas model 
within a statistical approach [32|]. The equilibrium statistical model was success fully 



used to describe the data on hadron multiplicities in relativistic heavy-ion collisions [33|. 
The system was treated in Grand Canonical (GCE) and Canonical (CE) Ensembles with 
taking into account the conservation of baryon, electric and strangeness charges. The 
decay of resonances was included in fiuctuations which were calculated for nuclear states 
to be realized along the experimental freeze-out curve. The whole available energy range. 
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from SIS to RHIC energies, was covered in this study 32|. 

We exemplify these predictions by the results in Fig{Tn] for the scaled variance of 
negatively and positively charged hadrons, 



--< (AiV")2 > / <N- > and uj+ =< {AN+)'^ > / < N+ > 



where A^^ is the number of particles in the given event and < ... > denote averaging over 
all events. 
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FIG. 16: The colliding energy dependence of the scaled variances for negatively uj~ (left panel) 
and positively lo~^ (right panel) charged particles calculated along the chemical freeze-out line 
for central Pb+Pb (Au+Au) collisions [32]. Different lines present primordial and final (i.e. 
including resonance decays) GCE and CE results. Solid points show the actual energy values of 
AGS (triangles), SPS (squires) and RHIC (dots) . 



As is seen from the comparison of the GCE and CE results in FiglTGl the exact charge 
conservation is very important in the Nuclotron energy. While the energy Eiab ~ lOA 
GeV is approached, the resonance decay is getting sizable. There is essentially differ- 
ent behavior of the scaled variances for negative, uj~ , and positive, uj~^, charges in the 



Nuclotron energy range Eiab ^ lOA GeV 



32|. 



At a first glimpse, these statistical results may serve as a reference point to look for 
peculiarities related to a possible formation of the mixed phase. However, it is not the 
case due to at least two reasons. First, in addition to the considered statistical fluctua- 
tions, heavy-ion collisions generate dynamical fluctuations. Secondly, the above presented 
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results correspond to an ideal situation when all final hadrons are accepted by the detec- 



tor. Indeed, as was shown in 



34| with using a transport approach to heavy-ion collisions, 



the fluctuations in the initial energy deposed in the statistical system yield dynamical 
fluctuations of all macroscopic parameters. In its turn, these fluctuations are dominated 
by a geometric variation of the impact parameter. However, even for the flxed impact 
parameter a number of participants fluctuates from event to event. The centrality of the 
selected events is commonly controlled by a number of measured projectile spectators. It 
was noted in 3J| that this procedure introduces some asymmetry between fluctuations in 
target and projectile, thereby the scaled variances uj~ and uj^ behave differently in the 
target and projectile rapidity regions. Thus, the above noted flrst and second reasons 
are turned out to be closely related. A flrst attempt to compare these results with the 
NA49 data at Eiab = 158A GeV shows a remarkable disagreement but unfortunately it 
is still impossible to make any conclusions from this fact due to large uncertainties in 
these experimental data 3J]. Note that the considered approaches take into account no 
phase transition. To search for signals of the mixed phase in fluctuations, the observable 
effect should be tested against its variations in bombarding energy, centrality selection 
and isospin asymmetry of colliding nuclei. 

One should stress that measurements of fluctuations require tracking detectors of large 
acceptance, good particle identiflcation and a precise control of collision centrality on 
event-by-event basis. Previous experiments suffered either small phase-space coverage or 
limited tracking and particle identification, therefore new measurements at the Nuclotron 
energy are of particular importance. From the experimental point of view the Nuclotron 
energy range seems to be ideal for these measurements. This is because moderate particle 
multiplicity and their relatively broad angular distribution simplify an efficient detection 
of all produced charged particles. 

A particular property of the mixed phase is the so-called distillation effect: While the 
total charge is conserved, its distribution between two phases is different. The structure 
of the mixed phase is especially rich and complicated if the conservation of the baryon 
number, electric charge and strangeness are taken into consideration simultaneously. The 



distillation effect is expected to result in certain observable effects as discussed in [35 



One should note that the use of neutron-rich isotopes or, generally, radioactive beams 
has recently become a central theme in nuclear and astro-physics researches. At present. 
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there are several facilities devoted to nuclear physics studies using low-energy beams 
of radioactive species including the Spiral at Ganil (Caen, France), ACCULINNA and 
DRIBs (Dubna, Russia), Cyclotron Research Center (Louvain-la-Neuve, Belgium), TRI- 
UMF (Vancouver, Canada), Holifield Radioactive Ion Beam Facility (Oak Ridge, USA). 
These capabilities for radioactive beams will eventually allow for detailed studies of the 
structure of nuclei on the path of astrophysical r-process and provide for fundamental 
nuclear structure studies of very neutron-rich nuclei but with small atomic numbers. 
Heavy-ion radioactive beams are also of great interest for both nuclear structure studies 
and some astro-physics problems, in particular, clarification of the role of new i^p-process 
in nucleosynthesis of nuclei with A ~ 100 37|. To our knowledge, there is no project to 
get accelerated heavy neutron-rich nuclei. Unfortunately, heavy nuclear isotopes which 
can be obtained in a reasonable amount cannot really be extended beyond the region of 
the isospin asymmetry parameter 0.39 < Z/A < 0.41 reachable with stable nuclei. 



IV. CONCLUSIONS 



A study of the phase diagram of strongly interacting QCD matter in the domain 
populated by heavy-ion collisions with the bombarding energy < lOA GeV and a search 
for manifestation of the mixed phase formation seem to be a very attractive task. The 
use of the isospin asymmetry as an additional conserving parameter to characterize the 
created hot and dense system draws new interest in this problem. Unfortunately, the 
available theoretical predictions are strongly model dependent giving rather dispersive 
results. There are no lattice QCD predictions for this highly nonpertubative region. So 
much theoretical work should be done and only future experiments may disentangle these 
models. 

A JINR Nuclotron possibility of accelerating heavy ions to the project energy of 5A 
GeV and increasing it up to lOA GeV can be realized in about two-three years. This will 
enable us to effort a unique opportunity for scanning heavy-ion interactions in energy, 
centrality and isospin asymmetry of the system to search for the mixed phase of QCD 
matter. For the latter point it seems to be optimal to have the gold and uranium beams 
in order to scan in isospin asymmetry in both central and semi-central collisions at not 
so high temperatures favorable for observation of the phase boundary reduction. The use 
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of strongly deformed U nuclei is quite promising to probe the orientation effect in heavy 
ion collisions. All this gives a chance to address experimentally many recent problems 
within the next several years before the FAIR GSI accelerator comes into operation. Being 
supplemented by scanning in the isospin asymmetry parameter, as discussed in the present 
paper, the proposed research program at the Nuclotron [3] may be considered also as a 
pilot study preparing for subsequent detailed investigations at SlS-100/300 [1] and as an 
integral part of the world scientific cooperation to study the energy dependence of hadron 
production properties in nuclear collisions. 
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